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In the present work we investigate a new statistical ensemble, which seems logical to be entitled 
the open one, for the case of a one-component system of ordinary particles. Its peculiarity is in 
complementing the consideration of a system with the inclusion of a certain surrounding area. The 
calculations indicate the necessity of taking into account the surface that delimits a given system 
even in the case when the latter is a part of a uniform medium and is not singled out one way or 
another. 

The surface tension coefficient behaves unlike two-phase systems in equilibrium and depends on 
two variables - pressure as well as temperature - and belongs to the boundary separating a hard 
solid from a fluid. As for the mathematical mechanism ensuring the fulfillment of thermodynamic 
relations, the emphasis is shifted from operating with series, like in the grand canonical ensemble, 
towards employing the recurrence relations of a new class of functions that incorporate Boltzmann 
and Ursell factors as their extreme cases and towards utilizing generating functions. 

The second topic of discussion that the present article deals with is the consideration of the surface 
tension and adsorption observed at the boundary of a solid body and a liquid or gas carried out on 
the basis of the analysis of the classical system found in a field of force of general type. The surface 
terms are calculated with the aid of field functions and the correlation functions of an unperturbed 
volume phase and behave somewhat vaguely; particularly, as a function of activity, they may start 
with a linear or quadratic term. 

PACS numbers: 05.20.Gg, 05.20.Jj, 68.08.-p, 68.43.De 



I. INTRODUCTION 

The Gibbs's grand canonical distribution is one of the 
cornerstones of statistical physics and its basic part leaves 
no place for doubt. Nevertheless, it seems worthwhile to 
check how well the relations of the grand canonical en- 
semble (GCE) meet certain criteria. One of those criteria 
is a specific kind of conformity that consists in a subsys- 
tem distribution (calculated on the basis of the GCE ) 
matching the initial one. 

Right from the beginning, the situation asks for such 
a test. Indeed, a well-known study of nucleation that 
can be traced back to J.W. Gibbs 1, p. 242] does exam- 
ine the surface of a nucleus. However the original grand 
canonical distribution describes fluctuations of the par- 
ticle number but at first sight does not seem to include 
any terms that represent surface. 

The incorporation of the terms describing the interac- 
tion over a surrounding surface into the statistical dis- 
tribution requires considering the environment that lies 
outside the limits of the system under analysis. This 
brings forward the idea of the open statistical ensemble 
(OSE). It can be further developed in the course of car- 
rying out the conformity check. 

Take note that in this case surface terms emerge when 
we describe a uniform medium that has no real surface 
but only a hypothetical one setting bounds to the sys- 
tem under consideration. This unexpected result can be 
better understood if we pay attention to the fact that 
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the point of interest for us in statistical distributions is 
the fluctuations occurring strictly inside the limits of a 
singled out volume. The states corresponding to these 
fluctuations already contain the delimiting surface which 
causes the abovementioned effect. Surface terms show 
the largest values for large fluctuations and are balanced 
out for mean values. It is also quite curious that a sur- 
face term shall be part of a standard GCE distribution, 
though with some deformations. 

With the treatment of the OSE distribution as a fluctu- 
ation probability, coherence between the thermodynamic 
and statistical approaches is established and the need for 
including surface terms into the distribution is justified. 

In an ordinary case of a two-phase one-component sys- 
tem the surface tension coefficient depends on one vari- 
able only - either pressure or temperature - due to the 
existence of the phase equilibrium curve. The situation 
under consideration extends beyond the reach of this con- 
dition and the surface tension coefficient starts depending 
on two variables. As it will be later shown it corresponds 
to the surface tension observed in the situation when a 
fluid interacts with a hard solid. 

Lately, the notion of the surface tension existing at the 
boundary between a solid body and a fluid has been on 
the agenda. For instance, it has been developing under 
the framework of density functional theory. The present 
article, however, suggests an alternative approach based 
on investigating a classical system situated in a field of 
force of general type which yields the combination of the 
field functions (solid body) and the correlation functions 
of an unperturbed fluid array. 

This approach produces the results that coincide with 
the direction initiated in [2j; however, it is significantly 
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simpler than the latter one allowing us to write the gen- 
eral form for the expansion in powers of activity. 

As a result, we've obtained the expression for £1- 
potential that gives us, in the limits of smooth fields, log- 
ically justified thermodynamic expressions for a nonuni- 
form system. For the fields that are microscopically in- 
homogeneous we obtain the surface terms corresponding, 
for the potential of a hard solid, to the coefficient that 
is part of the open ensemble partition function. For re- 
alistic potentials the volume term, surface tension and 
adsorption are not likely to be subject to strict separa- 
tion within the limits of a transition layer. 

Such treatment, as a matter of fact, deals with an in- 
soluble (nonvaporizing) solid body, adjacent to a liquid or 
gas. It seems only logical that in the setting of nonequi- 
librium the equilibrium curve does not have impact on 
the result. 



II. CAVITY FORMATION PROBABILITY 

This quantity is a vital parameter that defines the par- 
tition function of an ensemble. Thus, we will start from 
dealing with the first and most important term of this 
distribution and then we will switch to other ones. Let 
us begin with the canonical ensemble, move forward to 
the GCE and eventually come up with the definition for 
the OSE. 

Consider the probability of a fluctuation hole forma- 
tion inside a uniform statistical system, that is an area 
with the volume v that is devoid of any particles. In or- 
der to facilitate the calculations we will, in the majority 
of cases, have the integration volume set by the window 
functions tpi that are determined by the following rela- 
tion: 
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where v indicates the area volume while the integration 
will be considered to be performed over infinite space 
unless stated otherwise. 

We should also note that the algebra of window func- 
tions plays a significant role in the approach being de- 
scribed. Apart from that, as it will be further seen, the 
importance of these functions is high from the point of 
view of determining correlation functions. 



A. Canonical Ensemble 

The required probability in an isolated system shall 
evidently be defined as 
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exp(—(3U 1 N _ N )dr 1 ...dr N . 

(2) 



Here N is the number of particles in the system, 
[3 = 1/kgT, kg - Boltzmann constant, T - temperature, 
Uy N is the energy of particle interaction, and V is the 
system volume. C at p v c indicates that the probability 
is defined for the canonical ensemble - on the condition 
of finding N particles inside volume V . Zn is the config- 
uration integral: 
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By expanding the product in |(5J) and separating the 
integration over v, we will have: 
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where P[ k has the form: 
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(k) 

Py k defines the density of probability to find a given 
space arrangement of a set of specific particles 0, p. 181]. 
Considering distribution functions for a set of arbitrary 
particles 



N\ 
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(fc) 



Sc,x...k - (iv-fc)! xi - fe ' 
from ((3]) we shall have: 
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Pc,o = 1 
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9c!i...k dr i- dr k- (7) 



The analysis of a canonical ensemble causes some in- 
convenience due to distorted asymptotic forms of the 
functions k 0, p. 408] , which results in integrals be- 
coming non-local. Thus, a shift towards a GCE is re- 
quired. 



B. Grand Canonical Ensemble 

Let us average the Eq. Q over the fluctuations of the 
number of particles by applying to both of its sides the 
operation X)iv=o wnere 



Pn = 



z N Z 



N 



NIRx 



(8) 



- is the probability that the GCE shall have a definite 
number of particles ./V inside region V. Here z - is activ- 
ity: 
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where fi - is chemical potential, A = h/\/2nmkBT, h - 
Planck constant, m - particle mass and Sy - is the large 
partition function of a system with the volume V: 
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During the averaging we may change the order of sum- 
ming in the right-hand side ([7]) which yields 
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Here 



(fc) 

e'c.i. 



or in another way 
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x exp(-pU£+£ +k )dr k+ i...dr, 



fc+JV- 



ggj ^("0^) " is t ne ^" partical distribution function for 
a GCE. As stated above, these functions set the density 
of probability to find a certain configuration of arbitrary 
particles. For an ideal gas Qq\ k = g k , where g — N /V 
- is the density of particles. 

ip v cannot be regarded as a range of definition of func- 
tions in a general sense of the term, since, as we will fur- 
ther see, coordinates of free particles may extend beyond 
the boundaries that it sets. For the sake of conciseness, 
we will be calling it a range of assignment. For instance, 
in Eq. (fl3|) ijj v is the range of assignment of Qq \ k . 

In this case these functions make their entry through 
the series in the numerator (|13[) as well as in the denomi- 
nator - through the series (jlpp for Sy . Further on, it will 
be analyzed more in detail. 

Further investigation will reveal how significant the 
range of assignment ip v is, and for this reason we will 
be explicitly denoting it like we do with the belonging to 
the GCE-type by indicating it with G. 

Usually, the range of assignment ip v is not paid enough 
attention, thus being barely registered and weakly con- 
trolled, but according to the treatment that this article 
offers it is vital. 

Let us take the logarithm of the Eq. (fTT|) in order to 
localize the integrals. 



their range of assignment and type affiliation whenever 
it is necessary. 

Some relations will be true both for functions of the 
GCE-type and the OSE-type ones. In these cases, we 
will be leaving out the attributes of the GCE to avoid 
the exceeding sophistication of formulas. 

We shall now consider the properties of correlation 
functions and factors before moving any further. 



III. CORRELATION FUNCTIONS 

We shall further need the correlation functions of three 
types: complete, localized and partially localized ones as 
well as the three corresponding types of factors. 

Aside from that, the correlation functions will vary ac- 
cording to the type of the range of assignment. 



A. Complete Correlation Functions 

To draw more certainty to the definition this term will 
be applied to functions Qq\ j e ('^ V )j determined by Eq. 
(fT3]) and to their contiguous functions from the OSE - 
fc; which will be defined further. 



We will use their following property: 



QG ) ,i...k(^ V ) dr k 
d 



(15) 
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+ N-k+l 



This recurrence relation can be easily obtained from, 
for instance, Q12p by differentiation performed with re- 
spect to chemical potential Q. 

We shall also note the following property of their com- 
plete integrals: 



0G\...k('>P V ) dr l-~ dr k 



(N-k)) 



, (16) 



which follows from the definition. 

Equations (fTB"|) and (fTB)) are precise; however, the very 
values (N) for the GCE, as we shall witness further, in- 
corporate inaccuracies connected to surface terms. 

An important property of these functions is the split 
into a product as two groups of particles are drifting 
apart: 



til...k+l 



(fc) (0 
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(17) 
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.1=1 



^F { G ) 1 ... k (V)dr 1 ...dr k . 



(14) 

Here Fq\ fc('0 V ) - are Ursell functions or localized 
correlation functions We will be also pointing out 



B. Localized Correlation Functions 

Localized correlation functions or Ursell functions 
(they are also termed connected, inner and truncated cor- 
relation functions) have been analyzed by many authors 
and, particularly, in publications [6H9( . 
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These functions exhibit the local property for they 
quickly die away as any groups of particles, including 
a single particle only, start distancing from each other. 
This is observed due to property (|17[) . Let us write some 
of the first functions for a uniform system. 



F (4) - n [A) -(n {3) + n (3) + n {3) 4- n {3) )n 
1,2,3,4 — ^1,2,3,4 1^1,2,3 ' ^2,3,4 ^ ^1,3,4 ' Vl,2,4Jli 

For an arbitrary k > 1, we shall have || Chapter 9]: 
^ fc) fc = E(- 1 )^ 1 ^-l)!n^ ) ({"4) 



{»} 
1< k n < k 



a=l 



(19) 



where {n} stands for a certain decomposition of a given 
set of k particles located at r%, ...r k into non-overlapping 
groups {n a }, with I being the number of groups of par- 
ticular partition, k a - the size of group number a, with 
the summing performed over all possible divisions. 
An inverse relation of the following type is correct: 



{n} a=l 

The (fTTj) <-» ([14]) type of relation is universal. If 



(20) 
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^dr^.drk (21) 



and 



then 



K« = ln[u(«], 



fc=i j 



n« 



(22) 



(23) 



where AH t and w k \ are sets of functions bound by the 
relations of types JTl]), ((20j) 0j- 

For T^\ k {ilj v ) the following recurrence relation is 
true: 



d 



d/3fj, 



-k+1 



4*r!l-i(^)- 

(24) 



Obviously, (|24|) can be obtained from (|T5|) on the basis 
of definition HH). The fact that (N) vanishes from {24|) 
as opposed to Eq. (fT5|) reflects the local property of these 
functions indeed. Eq. (12"4")) is equivalent to 



(18) / ^^ G k l. k (^ V )dr k = z k §- z 
Multiple integration leads to: 
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^G\... k (' l P V ) dr l--- dr k 
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(26) 
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where I = 2, 3...fc. For I = 2 we shall write 

it 



V 



^GL k ^ V )dr2-dr k (27) 



dz 



fc-i 



Taking into consideration that 
'dP 



Q 



we eventually have 



= zp(—) 
dfi J T \ dz J j, ' 



IH 1 
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i...*(V' v )dr2...dr fc « z k (3 



Qkp 

dz k ' 



(28) 



(29) 



Eq. ([29]) , that initially appeared in Q during the inves- 
tigation of the analytical properties of P(z), is, generally 
speaking, approximate because qg does not correspond 
to q in the vicinity of the boundaries of volume. This 
question will be studied further on. 

We shall also note that integrals (|27p and (|29p . as a 
consequence of the local property of inner correlations, 
are defined by the microscopic region of integration in 
the vicinity of the first particle, unless the rank of the 
functions enters the macroscopic domain k ~ 10 23 . 

These are equations (l24l) - (|27)l . (|29|) that are analogous 
to Ornstein-Zernike relations for the GCE rather than the 
integrals found in |8]. 

C. Partially Localized Correlation Functions 

These functions constitute a hybrid of full and localized 
correlation functions. They play an important role in 
the mathematical apparatus of the OSE. As far as we 
are concerned, they were first introduced in work for 
the expansion of complete correlation functions. We will 
further see, though, that it is more convenient to employ 
the correlation functions normalized by z m+k , where m + 
k is their total rank. 
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Some of the particles that are part of these functions 
behave like those belonging to a complete correlation 
function, i.e. they do not cause dying out as they move 
away, whereas the other part shows local properties. 

Let us set partially localized correlation functions as: 



A 



(30) 



where superscripts m and fc set the number of delo- 
calized and localized particles correspondingly (to = 
1,2, 3, . . . , k = 0, 1,2, . . .). In their turn, the subscripts 
denote the coordinates of the particles with the restric- 
tion that the first m particles are treated as delocalized, 
and the rest of them as those possessing local properties. 

The structure of these functions is analogous to that 
of localized correlation functions of the (fc + l)-th rank, 
on the condition that the first to particles (the delocal- 
ized ones) are considered as one compound particle when 
writing (fTi?)) . 

In other words: 



(31) 



(j-i)in* 



(fc a +(m-l)5 a7 ) 



(K}) 



Q = l 



1< k a <k + l ^ k a = k + 1, 



where the denotations are analogous to those from (fT9|) 
on the condition that the summation is performed over all 
possible decompositions of fc + 1 particles, one of which 
is compound. S aj is the Kronecker delta, and 7 is the 
number of the group containing the compound particle. 

We may further write some of the first functions (over 
a local grouping) for a uniform system: 



4 m f =Q^.Jz m (32) 

J\i m ^) - \J m+1 ) _ J m ) n\/z m+1 

Am.2) _ r (m+2) _ (m+1) _ (m+1) 
• /H l...m+2 — Lfl...m+2 t'l...m+lt' fl...m,m+2^ 

t'l...mt'm+l,m+2+ z t'l...mP J/ z 



The first functions over the delocalized group coincide 
with the completely localized ones with the accuracy up 
to a multiplier: 



A 



(i,fe) 

l...k+X 



T (k+1) , k+1 
• r l...k+l' /i 



(33) 



Equations (fT!))) , (j3"Tj) - (f3"3")) show us that, evidently, 

functions A ™ '^ +k , in a certain sense, fill all the space 
between complete and localized correlation functions. 

The counterpart of ([25]) can be easily obtained from 
the definition (f3"Tj) : 



^+k^ [ G,ll m +k^ V ) dr m+k 



u ^GS...m+k-l\V ) 



dz 



Hence follows the counterpart of (f2l)]) : 



(34) 



n ^ 

. i=! 



where I = to + 1, to + 2, .... to + fc, and for (|27| we write 
the analogous expression: 



" 4 G?i fc) m+fc(^ y ) dr; -- dr ™+ fc ( 35 ) 



n w 



_z__m+l 



(36) 



or, taking (|32p into consideration yields in: 

m+k 



n ^ 



_2__m+l 



^G?i fe) r™+/c(V' y ) dr m+l--- dT "m+'c (37) 



gk 



We shall also underline that in order for these rela- 
tions to be true integration should be only carried over 
the coordinates of localized particles. The fact that the 
degrees of activity vanish from (f3"4")) - (f3"6"]) in contrast to 
(|25| - (|27| gives us an insight into the purpose of the 
normalization. 

The generating function for \ is the logarithm, as 
we can see from pT|) . (fl~4"|) and ([22")) . The generating 

function for -4x m+& " * s * ne fraction of the following 
type: 
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. i=l 



(38) 



4 m .m+/c dr l-- rfr ™+fc- 



As it can be easily obtained from (|38j) by multiply- 
ing the series this equation is ensured by the recurrence 
relation 



Am.k) Am+k,0) 
■^l...m+k ~ ^l...m+t 



(39) 



EV /|( n > ) A(m.,k-n) 
/ j •^ t l...n vH n+l...m+)fe' 



n—1 perm 
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where we employed the expression for A^ n ' ^ from (|32f . 
while the internal summing is performed over the permu- 
tations of localized particles only. 

Relation (|3"9")) can be proved either by means of direct 
exhaustion of all decompositions in accordance with <|3 1 1) , 
or by the way of a repeated substitution of the expression 
for A^ n '^ J _ k in the right-hand side of (|39p. 

Work [7j features another recurrence relation in regard 
to Ursell functions that, upon translating it to the lan- 
guage of ^4^ m '^ +fc , may be written as: 



• /H i...fc+i 



•^l.-.fc+i 



fc 

£ 

71=1 



(40) 



Vl+1. 



n) 

.fe+1 



where the brackets stand for the averaging over the per- 
mutations of particles. 

As a matter of fact, it is the very same equation (j3"!)j) 
at to = 1 but written in a symmetric form. The question 
about the necessity of the symmetrization of all recur- 
rence relations for Ai'^, k , which we will further face, 
falls beyond the scope of the present work; note, though, 
that in the given case there is no need in symmetrization 
as Eq. ([55]) is strictly satisfied in the assymetrical form. 
Nevertheless, (l3l?l) can also be written in the symmetric 
form, but this will require the averaging of the left-hand 
side of the equation and not that over the sum permuta- 
tions only. 

The satisfaction of these and many other recurrence 
relations for •A^ n '^ +k is crucial for the mathematical sup- 
port of the approach being presented. As it will become 
clear further in the course of our discussion, they are also 
true for the corresponding factors and give us the oppor- 
tunity to operate with the expressions from the OSE. 



D. 



Factors 



Replacing g\ ' fe by z k exp(— j3U* k ) in the defining 
equations p^|) . (f3"Tj) brings for the following factors: 



Boltzmann's: 



Ursell's: 



exp(-/3E#„ fc ) «-> g[ k) . k /z\ 



and partial localization factors: 



j(m,fe) j\(m,k) 



„yrn, K) ,yn 
a l...m+k "^...m+fc- 



At the same time 



1 <-* g/z. 

In other words, let us define B 



(m,fc) 
l...m+fc 



(41) 



(42) 



(43) 



(44) 

by an equation: 



(45) 



x nexp(-/3C/ fc =+( m - 1 ' 5 -({n Q })) 

a=l 

1 < k a < k + 1 , a =/= 7 or m = 1 

I 

1 < fc 7 < k + 1, m > 1; kg < k + 1, 

Q=l 



whose notations are analogous to those of (f3"Tj) . The dif- 
ferences are observed due to the fact that one-particle 
groups do not have any contribution to the products in 
this case. 

As factors are primary values independent of the en- 
semble type or window functions, we omit these notations 
for them. 

The generating function for B^'^+k' as follows from 
(|3"5)) . (fT3"|) , (|4"Tj) is a complete correlation function of the 
GCE-type 



(m) 
6(3,1.. 



£ AT J 
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n ^ 

_i=m-\-l 


e l3u ™^+ k dr m+1 ...dr m+k 


00 k p 

k=l J 


' k 

IK 

.7=1 


e pu *- k dri...drk 



(777,0) 



-T- 

k=i 



n # 



,1=771+ 1 



(46) 



Ursell factors, also known as cluster functions, are in- 
volved in a well-known expansion of pressure in powers 
of activity 0, p.129], 0, p.232]. 

Boltzmann factors are also well known and make part 
of configuration integrals ^ . 



The factors of partial localization, however, have not 
yet been investigated - at least, as far as we are con- 
cerned. 
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For instance: 



UP 



u 



= 1 

= cxp(- / 3C/ 1 2 2 
t 

r 2 



(47) 



i (3 L - exp(-pUl 2j3 ) - exp(-/?E7? i2 ) 



exp(-/3Etf 3 ) - expi-PUia) + 2 



The first Bx!?.'m+fc' s over ^ ne localized group: 

= exp(-/3f/ 1 m . ro ) (48) 
= exp(-/3£/££ +1 ) - exp(-/3C/™ J 

= ex P (-/3C/r. + ™ +2 ) " exp(-/3f/™^ +1 ) 
- exp(-/3C/™+^ /m+2 ) - exp(-/3£7™ J 
x exp(-/3tf» _+.,) + 2exp(-/3E/£.J 



and over the dclocalizcd one: 



(i,fe-i) 



(fc) 



including 



(1,0) 



1. 



(49) 



(50) 



where we showed the form for a uniform environment. 

As it will be later shown, the factors constitute the 
limits of the corresponding correlation functions that are 
normalized by the required degree z with activity tending 

to zero. In their turn, A.^ m+fc a ~° n0 ^ rec i u i re normal- 
ization in order for this to be true. 

In the same manner to A^ n '.^ +k , factors B[ m '^ +k gen- 
eralize Ursell and Boltzmann factors incorporating them 
as extreme cases. 

Evidently, factors maintain the algebraic properties of 
correlation functions and thus the recurrence relations 
identical to (l39l and (l40l) hold true for them as well: 



(51) 



n+l...m+fc' 



„{m,k) 
°l...m+fc 



r>(m+k,0) 
°l...m+fc 
fc 

£j(n,Q) j g(m,fe-n) 

n— 1 perm 



true when to > 1, k > 1 and 



pj(l,m) _ „(m+l,0) 
°l...m+l — D l...m+1 



(52) 



E 



B(n,0) K?(l,?n — 7i ) 
1_ ^ &n 



n+l...m+l 



a number of cases we will make use of the symbol s = ; at 
that we will imply that the right-hand and the left-hand 
sides of the equation underwent the symmetrization over 
those indexes where it is necessary. In this case, sym- 
metrization will not only mean the averaging over the 
permutations but also aligning the indexes to the natu- 
ral scale if it is required. This operation appears as a 
consequence to property (|50[) . 

So, apart from (|51[) . ([52")) the following equations are 
fulfilled: 



.k) sim 

7n-\-k 



*j(m— 1,/c+l) 



(53) 



fc 

E 

71=0 



B 



(m — l,fe— n) *2(l,n) 



l.,.m+fc- 



rv ,B 

n— 1 ■ 



m+fc— n. . .m+fc ' 



at fc > and m > 2, as well as 
,(i,fc) ^ /'fc 



B 



.fc+i 



fc 

E 

m — 1 



X 



(-i) r 



(54) 



j2(m+l,fc— m) 
e l...fc+l 



perm 



where fc > 1; and we also witness the satisfaction of two 
classes of recurrence relations that have two free indexes 
with the former class containing the low-value indexes of 
delocalized particles: 



k 

E( 

m=0 



i \m ( '" \ yj(m+l,t+k — m) sim 



E 



(55) 



(-D'E 



n=0 



K?(/,fc+n— /) — n) 

"l...fc+n "fc+n+l...fc+t+l' 



where fc > 1, t > 0, and for fc = 1 it reduces to (|53"|) with 
rn = 2, while for t = it amounts to (|54[) . The latter 
class, au contraire, contains those of localized particles: 



EM) 

=0 

x(-l) 



m=0 



„(m+t+l,k—m) sim 
"l...t+fc+l 



t 

<E 

n=l 



1=0 



r>{l-\-n.k — l) — n) 

°l...fc+n °fc+ri+l...fc+t+l 



(56) 



(s,k-s) H (l,t) 

fc+l...fc+t+l' 



where fc > 1, t > 1, and for fc = 1 it reduces to ([52)1 . 
Equations (f53"|) - ([56")) will be proved in the section fV Dl 
All of them ensure the fulfillment of certain physical 

relations that will be defined later. 



where to > 1. 

As it has been already said, the question about the ful- 
fillment of recurrence relations in the asymmetrical form 
of type (f5Tj) extends beyond the frame of this work. It 
is not improbable that all of the further given relations 
have this form but we will be using the symmetric one. In 



IV. THE FORMULATION OF THE OSE 

A. The Partition Function of the OSE 

Now we are able to compute pjj and consequently define 
the partition function of the ensemble. This function 
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will further be referred to as Y„, where v stands for the 
volume of the system: 



t„ = 



1 

p5" 



(57) 



Let us employ the expansion of functions !Fq\ k(il> v ) 
that is given in Writing it in the language of ip v 

gives: 



^L fc (tf v ) = *^£!*+**E 



(58) 



k-\-n 

n $ 

l=k+l 



By substituting (|58f in (|14[) and changing the order of 
summing we have: 



(59) 



t-i ' ^ 

\k=l ^ ' i=lj=k+l / 

Recasting (|5^|) with the aid of the binomial formula, 
and utilizing the symmetry of U^f 1 t with respect to par- 
ticle permutations, results in the expression: 

°° t t- 



II 



Ul' t dr 1 ...dr t (60) 



£5 



Ui) t dr\...dr u 



or, 



(61) 



Indeed, both terms in the right-hand side of (|6"0|) con- 
stitute the logarithms of the partition functions of the 
GCE, since, taking the logarithm of Eq. (fTU|) . and tak- 
ing ([3]) into account as well as relation ([22]) yields: 



t=l 



hlf} t dri...dr t 



(62) 



Volume terms of Eq. (|62p give us the expansion for 
pressure [H, p. 135] but we will extract them a little bit 
later. The formulation of the OSE is concerned with the 
following limit: 

mY„ = lim (InSy — lnSy_ t) ) , (63) 

while the intensive parameters of the environment and 
the parameters of volume v are maintained, or course. In 
order to proceed to the limit we unite the integrals from 



the right-hand side of (I60|) and assume that ipf = 1j f° r 
all i: 



00 r 



i-ri( 1 -^) 



t dri...dr t . 



(64) 

Thus, we obtained one of the expressions for the par- 
tition function of the OSE in a uniform environment: 



{°° _t r r * 
Efr/ 1 n 
4=1 ' J L i=l 



-TTa-V'n 



Ui t drx...dr t 



(65) 

Now let us split it into volume and surface parts by 
reconstituting (|64|) in the following way: 



00 z t r 



+ E 



f! 



1 - - nt 1 - 



(66) 



Ui}_ t dr\...drt- 



The first term, being the volume one, brings forth the 
well-known expansion for pressure: 

°° t r 

Pv = fcsT^ Z - J il>lU^. t dr 1 ...dr t (67) 



or: 



P{z,T) = zk B T + k B TY j ^ j U [ ^ k dr 2 ...dr k , (68) 

fc=2 ' ^ 

it has, however, the counterpart for surface tension writ- 
ten below. Let us factor out 1 — ipl from the second term 
and replace it by Xii where 



x ■ = i - C 



(n e v) 

1 (ri i v). 



(69) 



If function localizes the i-th particle inside volume 
v, then Xi does the same outside its limits. We shall 
further have: 



lnT„ = Pv/k B T 

OO + 

t=2 



i-ri( 1 -^) 



(70) 
Ui} t dr\...drt- 



The second term in the right-hand side of ([70|) is pro- 
portional to the surface area. Indeed, its structure pre- 
scribes at least two particles to lie at the different sides 
of the boundary. If we expand the product in ([701) , then 
all terms will be written: 



i=2 



Uf' t dri...dr u 



(71) 



where 2 < j < t. When the first particle is fixed, the 
integrals of type (171]) are defined by the local area in the 
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vicinity of this particle due to the local character of Ursell 
factors. As a consequence, J ...dr 2 ---dr t do not depend 
on displacements of the first particle along the boundary 
of the system. However, when the first particle moves 
away from the boundary they quickly die away because 
of the fixing factors ip? and the local character of Ursell 
factors. Performing the integration alongside the surface 
and taking the area out of the summation sign yields: 



T„ = exp/%P(z, T) + aa(z, T)], 



(72) 



where a - is the area delimiting the system, vP(z, T) is 
given by Eq. (|67|) . aa(z, T) shall be written 



aa(z,T) = fcflTy^ 



4=2 



(73) 



Xi 



l-n^-c 



i=2 



Ui ' t dr\...drt 



and for 17(2;, T) we have the following expression: 



<r(z,T) = k B T^2 — 

t=2 



(74) 



Xi 



ri( 1 -^) 



Ul \ t dxidr 2 -.-dr t , 



where x\ - is the coordinate perpendicular to the surface 
given by the boundary of ip?. The direction of axis x\ is 
chosen in the way to satisfy dx\ > 0. 

The first terms of expression (174")) were obtained in 
by means of the diagram technique, but building up large 
number terms in the series with the aid of this approach 
seems to be extremely difficult. 

For the cavity formation probability we shall write: 



pi = cxp -0[vP(z, T) + aa(z, T)}. 



(75) 



We will discuss this result from the thermodynamic 
point of view later on. 

By making use of x"i t ne expression (|64|) takes the 
following form: 



00 1 

^ = Efr 

*=i 



1=1 



u[ t \dri...dr t 



(76) 



which is interesting to compare with the expression for 
the logarithm of the partition function of the GCE (p2"|) . 



of the zero term given above. 

V 

PC,m 







mj 


Zn 


n l 


n 



.N) 



(77) 



n n 

i—l j—rn+1 



dr\...dr 



N ■ 



The binomial coefficient appears in (|77j) as a conse- 
quence of us considering only specific sets of particles at 
this stage. By calculating in the same manner to (13}, ([7]) 
we have: 



N-n 

— y 



(-1)* 



TO! a — ' fc! 

fe=0 

m-\-k 



(78) 



(m+k) j , 
eC,l... m +k dr l- dr ™ + k- 



Shifting to the GCE and changing the order of sum- 
ming yields the expression analogous to (JTTJ) : 



PG,m — 



E 



(-1)* 



TO! * — ' fc! 

k=0 

m+k 



(79) 



Sai. k) rn+k(V)dr 1 ...dr m+k . 



We further divide the series (|T9[) by (|TT|) and thus ob- 
tain the relation: 



■ ! Pg,t 



E 

fe=0 

/ 



(-*)* 
fc! 

m-\-k 
. i=l 



(80) 



■^Gxlm+k v )dr!.. .dr m +k , 



where to = 1,2, ... We resorted here to the fractionary 
generating function (f3"B"|) for the GCE-type correlation 
functions. 

Now we shall derive the expression analogous to (I58[) . 

which generalizes it for ■A^Q^\ n+k (ijj v ). 

In order to do this we shall expand this function in a 
Taylor series in powers of activity 



(m,fc) 



V — 



Q z n^G,l...m+k\ Z ' V ) 



(81) 



z=0 



B. The Distribution of the OSE 

We can now calculate the common term of the distri- 
bution of the OSE. So, consider the probability to find to 
particles in volume v for a uniform system - j> v m . The out- 
line of operations here will correspond to the reckoning 



and taking (I35| into account yields: 

■ /X GS...m+k\ Z i V I " • /t G,l...m+)tl U ' V I 

m+k+n 



E 



n # 

i—m-\-k-\-l 



xA al k M+n(°^ V ) dr rn+k+l---dr 



(82) 
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Then, the definition (1131) , at z ^ 0, leads to 

^ fe) fe ^z fe exp(-/3C/ 1 fc .. fc ), (83) 

[3, p. 410]. As decompositions (f3Tj) always maintain the 
total rank Qi k we may conclude that 

+ *(°) - BjKU ( 84 ) 
and thus we obtain the relation we initially sought for: 



- A G,l...m+fcV 2 > V J — "l...m+fc + / , „| 

n v, 



(85) 



S(m,k+n) , , 
l... m+ fe +I j« r m+fc+l---a^m+fe+n- 



At A; > 0, m = 1, when taking (J33J) and (J49J) into 
consideration, we obtain ((58]) : at fc = 0, to > 1 and 
by taking ([3"!?]) into account it yields the expansion for 

&..fc(* V ), 



(86) 



n # 



,2— m+l 



n=l 



l...m+n m+l---^' m+ni 



which exactly agrees with (|4"6")l . 

Finally, at fc = 0, m = 1 the well-known expansion for 
the density of a uniform system follows from (185|) if we 
take (l49l) and (l50l) into consideration: 



n+l 
,i=2 



(87) 

(Evidently, (J57J) may also follow from ((SSJ) at fc = 1). In 
this relation the factors ^ are often neglected or gG(ip V ) 
is substituted for g, or both. 

Now we can substitute ([55)1 in ([SO]) , and, by following 
the steps similar to proceeding from (TTJJ via (|5§|) to ([GlT]) 
we obtain: 



Pb,r 



E 



" t=o 



m m+t 



(88) 



n n *x®Y-+i) 

i—1 j — m+1 



where to = 1,2.... The passage to the limit 

V^oo, (89) 
u = const, 
ip = const, 



where ip are the intensive parameters of the environment, 
yields 



to!T„ ^ t 

v t=o 



T- 



(90) 



m m+i 

n n « 



B ( i m :!rl +t dr 1 ...dr m+t . 



This expression solves the posed problem by defining 
the general term of the distribution of the OSE. 

It is easy to see that the first term of the series (|9"0")) 
corresponds to the distribution of the GCE ([8]) due to 
(|48p or because of ([50")) at m = 1 with the accuracy up to 
normalizing factors - partition functions. 

The structure of the series (|90l) seems to be of inter- 
est. The location of the delocalized group of to particles 
is set by window functions ipf inside the volume of the 
system. The location of the localized group of t particles 
is, in its turn, defined outside it by functions Xji being, 
however, closely bound with the volume due to functions 
B[ m 'm+f The summing is performed over all clusters that 
are bigger than to. It may be stated, that all volume- 
related properties of the expansion are assigned by the 
first to particles - ordinary, delocalized ones, whereas sur- 
face properties are predetermined by running t particles, 
that are localized. 

As the differences between the GCE and the OSE are 
related to surface terms, all the elements of the series 
([9H)l . starting from the second one, contain the volume 
as well as the surface of the system in different degrees. 



V. CERTAIN PROPERTIES OF THE OSE 

Making calculations with the OSE is somewhat more 
difficult than with the GCE. Every single operation with 
this distribution is ensured by a certain class of recur- 
rence relations for B^ m '^ +k . For brevity sake, we say that 
an operation yields a recurrence relation or class. On the 
other hand, the arising difficulties can be overcome in 
a number of cases by means of using various generating 
functions. Consider some of these operations. 



A. The Condition of Normalization 

J2Pm over au m ' s shall be equal to 1, and thus we 
obtain the equation: 



^— ' to! ^— ' 

m=l t=0 



(91) 



m m+£ 



n n « 

i— 1 j— m+1 



This is yet another form for the partition function of 
the OSE. It is difficult to see here expression (|65|) ! In 
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order to prove these formulas identical we may take their 
logarithmic derivatives. Then we have: 



1 dT v df3{Pv + aa) 



T v dz 



dz 



(92) 



B. Mean Particle Number 

The next operation to be considered is calculating to. 
It yields the class of recurrence relations ([5"5j) . We will 
proceed from the natural assumption that 



dTv =T d(3(Pv + aa) 
dz ' dz 



(93) 



moreover, T„ in shall be defined by expression (f9"Tj) . 
while Pv and aa - by (|67|) and ([73]) correspondingly, 
where t are substituted by 0\ ' 1 ^ according to (|49l : 

OC j - 

«P(z,T) = k B Tj2j i> v iBt* t - 1) dr 1 ...dr t (94) 



aa(z,T) = k B T^2 



(95) 



t=2 



Xl 



l-na-v?) 



B{ 1 : t - 1) dr 1 ...dr t . 



Let us change the order of summing in the expressions 
to facilitate the calculations. By the example of (|9"Tj) we 
write: 



T 



1+ Et7 / E 



(96) 



n n w 

i—l j — m+1 



Differentiating and multiplying the series as well as 
equating the expressions at different degrees of z and 
products ip? yields recurrence relations (fB"3")) . (|54"|) and 
the class of relations (|56|) . 

By applying ordinary thermodynamic formulas to (|92|) 
we may write: 



d In T„ 
d In z 



TO h - m s , 



(97) 



where 



TOfc = v 



dP 



is the number of bulk particles, whereas 

da 



to, = —a 



(99) 



- is that of the surface ones. 

We will further compare (|97|) with a similar relation 
for the GCE. 



to = nib. 



(100) 



The proof can be given by taking the following steps. 
According to (|90|) to takes the form: 



1 z m 

m = E r m _ i)i E 

m=l y ' t=0 



oo f 



(101) 



n n « 

i— 1 j— rn+1 



M m .m+t dT "l-- rfr m+*- 



Let us change the order of summing: 

m = y; E (7tt)T / E ( m _ i 

s— 1 rn — 1 v 



(102) 



n n * 

?—l j=m+l 



"X? 



6','";; '"V/r,...'/r- 



We shall equate this expression to (|98|) . where we need 
to employ (|94p for P. By multiplying the series and 
equating the expressions at different degrees of z and 
products ij>V we obtain the class of recurrence relations 

(USD- 

Equation (| 1 1 [) yields a new expression for m s . With 
substituting (|100[) in (|9T|) . we have: 



9 In T„ 

— = TO - TO S 

omz 



(103) 



Further on, by using T„ in the form of (|91[1 and to in 
the form of (|101[) we can obtain from (|103|) the following 
expressions: 



OO m OO 4 



n n w 

z— 1 j—m+1 



' m+t, 



or 



1 °° z s /" s_1 

^^J^Ty.J E 

s— 2 v 7 m — 1 



s- 1 

TO 



(105) 



n n * 

z— 1 j—m+1 



i Aj 



i; s 'dr 1 ...dr s 
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Correlation Functions of the OSE 



As we shall further witness, the OSE may utilize cor- 
relation functions of different types, including those of 
the GCE-type, but this section will deal with the corre- 
lation functions setting the probability of a given config- 
uration of particles in an ensemble, which are analogous 
to (fl3j). and which will be denoted, as it has been previ- 
ously stated, as Q{\- We will also consider here J-^ k \ 

and Aj^'^ +k developed on their basis. 

Let us turn back to Eq. (|56")l . As it becomes clear from 
(|63)) . (f90|) . the essence of the shift to the OSE is in the 
passage to the limit (|89|) . So, we may set 



lim g l ™l 



(106) 



thus, as a counterpart to (|86j) we have 



(107) 



'5Z~T / B{ m .Z ) +n dr m+ i...dr m+n 



(let us remind the reader that when the window functions 
are not indicated, integration is performed in the infinite 
limits). Instead of ([55)) we write 



•M...fc - z u \...k 



(108) 



Y.ij U [ ^l n dr k+1 ...dr k+n 



n=l 



while instead of (l85l we have 



A (m,k) , s _ „(m,fe) 
m+kV* ) ~ u l...m+k 



(109) 



\ " 



~ZJ I t} l... m +k+n ar m+k+l- 



.dr. 



m+k+n • 



Obviously, quantities •7 r {' c ' ) fe and -4x m+fc comply with 



relations (| 19[) and (|31j) correspondingly, and the passage 
to the limit can be performed always, because the terms 
of all series are integrated over the coordinates of local- 
ized particles only. 

Equations (|107|) - (|109p convincingly show that unlike 



their GCE counterparts quantities Qi_ m , ^i K \ an d 

•^i^'m+k P ossess the property of translation invariance 
and by that reflect the properties of a uniform system. 

In particular, the difference between these pairs is sig- 
nificant in the cases where a group of particles 1, ...k is 
located at a microscopic distance from the boundaries of 
the volume set by function ip v . However, there is an- 
other option apart from this case, where free particles 
are found outside the limits of this volume. 

The expressions analogous to (|34|) . (f35|) . ([37]) shall be 
written as 



-(fe) 



A {m,k) dr 

■^l...m+k u ' r, 



tA/H l...m+fc- 



(110) 



I A m :"lkdr,..dr. 



Qm+k-l+lj^™, 1 -™- 1 ) 
m+k = a. ' (HI) 



Q z m+k-l + l 

where I = m + 1, m + 2, m + fc, and 



Am,k) , , 

■ /i l...m+k aVm + 1 --- arm + k ~ ~Q^k 



h 



(m) 
Ql...n. 



(112) 



Once again, we remind that the integration should be 
performed over the coordinates of localized particles only. 

Eq. (|110p is easy to explicitly obtain from (|109[) . or 
from by applying the passage to the limit (|89"j) . 

Equations (jllip and (I112[) obviously follow from (|110p . 

Atm = l (fTT0|) - (fTT2|) yield the OSE counterparts of 

U25J, © and (231 for ^ fc : 



• r i...fc u/ fe 



(fe-i) 
i...fc-i 

yk-l 



Qk-l+1 

T { ^ k dn...dr k = z k Qzk _ l+1 



T 



(l-X) 

1...1-1 

rl-l 



where I — 2, 3...fc. For / = 2 we shall have 



F < t k \.dr 2 ...dr k = z f 



(113) 



(114) 



(115) 



Again taking 



into account we eventually obtain 

Qkp 



3 r i k l k dr 2 ...dr k 



z K l3 



dz k 



(116) 



Eq. (|116[) unlike (|29|) is precise, since from (|10T[) it 
follows that q corresponds to (|6"8")) and (|2"8)) . Thus, we 
have obtained the exact analog to the Ornstein-Zernike 
equation in the form of relations (| 1 1 5[) , (|116p . 

Clearly enough, a direct analogue of (fTS")) for the OSE 
is impossible. Indeed, as the environment is limitless, 
equations cannot contain extensive quantities. The ar- 
tificial singling out of a certain scope of integration will 
lead to a forced limitation of the configuration of free 
particles. Leaping ahead, we shall say that this limita- 
tion is eliminated when integration is performed over the 
coordinates of all particles, like in (|126p . 

The adjusted counterpart of (fTS"!) will be the form (|110[) 
at k = 1. Indeed, by developing it in accordance with 

(EH), 



r (m+1) (m) i j 

[£»!.. .m+1 - Qi...mQ\ dr m+k 



<9z 



„( m ) 



(117) 



we witness that, in the given sense, it corresponds to (TH 
By substituting series (|107[) in (|1 17[) we can clearly see 
that this expression yields the recurrence relation (|53p . 
which, in its turn, will be proved in the following section. 
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D. The Use of Generating Functions 



Reconstitute (|90| to a slightly different form by sum- 
ming the series over \ v ■ We obtain: 



1 



i!T, 



Q i bl.. m U v )dri...dr m , (118) 



where m (x v ) is defined by expression (f8"6"|) . 

Here we face the case where all the coordinates of parti- 
cles 1, ...to that are of interest to us lie outside the bound- 
ary of the range of assignment of the correlation function 
set by functions \ v ■ 

First of all, note the similarity of expression Q118P with 
the common term of the distribution of the GCE © . In 
the following approximation 



(m) 



z m exp(-/3f7 1 m 



(119) 



they correspond to each other with the accuracy to a 
normalizing factor, that is, in the limit of low densities 
when the first term plays the major role in the expansion 
of a correlation function. (This means, in particular, that 
for low densities the significance of the surface tension is 
reduced.) 

The condition of normalization yields the expression 



oo 



m— 1 



1 



QG*L.m(x v )dri...dr r , 



By taking its logarithm we obtain 



(120) 



' ml 

m=l 



T. 



(m) 

G.l— m 



(X v )dr 1 ...dr. n 



(121) 

By making use of ([55)) and changing again the order of 
summing as well as regrouping the products of the win- 
dow functions we obtain Eq. (|76[) . Since we started from 
(|9"0")) . it means that we proved the recurrence relations 
(|53| , ([54]) and the class of relations (156)) . 

Then, by utilizing (|118p we average the expression 
to(to — 1)(to — 2)...(m — k + 1): 



(to — k)\ 



1 



~k-\-n 



l 



(122) 



e < GA^k+n(x v )dr 1 ...dr k+n . 



Let us divide the series in (|122[) by employing T t , as 
series (|120l) . Likewise to ([50)1 we have: 



(to — k)\ 



n— 



(123) 



.4 



G,l...k+ 



n {x v )dr 1 ...dr k+r , 



By substituting the expression with the window func- 
tions x v which is analogous to (|55|) in (|123[) we obtain, via 
repeating the already habitual procedure of changing the 
order of summing and regrouping the window functions: 



■y- 

t=0 



(124) 

Setting k = 1 in (|124j) yields the expression for the 
average particle number to: 



rn 



oo t . 

EtF r i B[ 1 % 1 dr 1 ...dr t+1 , 



(125) 



which, evidently, corresponds to (J94J), (|98|) . Thus we 
have found an alternative proof for Eq. (|100[) and, conse- 
quently, we proved the class of recurrence relations (|55p . 

The summing of series (|124l) in accordance with f|10T[) , 
gives: 



(to — fc)! 



(126) 



which corresponds to (|16p . Thereby, this equation is pre- 
cise for the correlation functions of the OSE as well as 
for the GCE. However, as we will further see, the values 
of average quantities in the OSE are not distorted unlike 
in the GCE. 



VI. A SYSTEM IN THE FIELD 

The present article will be dealing with the issue of 
a system of particles in the field only to the extent, to 
which it is related to the major topic of discussion that 
is the formulation of the OSE. 

Let us now consider a statistical system located in a 
field where the scope of the field is far smaller than that 
of the system and with the field quite far from its bound- 
aries. For such a situation it will sussici to study only 
the GCE as we will further see. In other words, we will 
be interested in what happens inside the system and not 
on its external boundaries. 

Furthermore, the foregoing discussion will be based 
upon the presupposition that function P(z) is analyti- 
cal and the corresponding series converge. This issue has 
been raised in a large number of works, e.g. (Tl| . So, 
the present work does not discuss this question and as- 
sumes all functions to have analytical character and the 
corresponding series to be in the state of convergence. 

Apart from that, we shall suggest that the poten- 
tial of the fields under consideration decays/rises quickly 
enough for the corresponding integrals to converge in- 
ward or outward the body. 
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A. Basic Equation 

Consider the configuration integral of a nonuniform 
system: 

r N 

Zn = J exp(-/?^ Ul - (3Ut. N )dr 1 ...dr N . (127) 



Here Ui is the energy of interaction with the field of 
the i-th particle. 

Let us introduce the functions 



def 

<Pi = Lp{ri) = 1 - exp(-0Ui), 



(128) 



that depend on the coordinates of the i-th particle, with 
their advantage being the ability to localize it inside the 
scope of the field. The further elimination of cxp(— (3ui) 
from (fT27l) yields 



Let us take the logarithm of Eq. (|134|) . 

(-1)' 



InSy - InSy = ^2 ' 



fc=i 



kl 



(135) 



x / (f l ...ip k T ( ^\ k {ip v )dri...dr k . 



As it was mentioned before, we are considering the case 
where the distance from the field scope to the boundaries 
of the system is large. On these conditions, the compar- 
ison of (|58| with (|108[) makes evident, as a consequence 
of the local character of Ursell factors, with a very fine 
accuracy that 



(136) 



that is why in the frame of these assumptions the follow- 
ing is true: 



Z U N = / (1 - ip x )...{\ - tp N ) exp(-puH. N )dri...dr N . 

(129) 

For an analogous system devoid of the field we will 
obviously obtain: 



Z N = J exp(-{3Ui N )dri...dr N . 



(130) 



lnS^-lnSy 



fc=i 



(137) 

This equation is basic for the further consideration 
of the field related issues. As we can see, the vari- 
ables here are also separated: ip are the functions of the 
field, whereas Ursell functions belong to a uniform system 
without a field. 



As for the fraction of the partition functions, by in- 
tegrating over the coordinates of free particles we may 
write 



7 v 

N 

Zn 



1 + E (I) ( ~ 1)fc / ^■■■9kP[ k) . k dr 1 ...dr k . 

(131) 

Moving on to the distribution functions for an arbi- 
trary set of particles Qq\ k from p31[) follows: 

7 U N (— 11* r 

^ = 1 + EHr / r:-^!>r A ^lr l ...dr,. (132) 



Let us turn to the GCE and average Eq. (I132j) over the 
fluctuation of the particle number for a system without 
a field. By introducing the value 



E 

7V=0 



A! 



(133) 



which is, evidently, a large partition function for a sys- 
tem in the field, and by changing the order of summing 
in the right-hand side we shall have 



. i - - i> « — 



i+E 



' (pi...ip k g^\ k (ip v )dri...dr k . 



(134) 



B. Smooth Field 

Let us set the scale where a field is subject to significant 
changes as macroscopic. Then, due to the microscopic 
region of the decay of we may employ f|l !€>[) and 

write 



1=2 



where k = 2,3.... (It is necessary to mention again that 
(j!38[) is true for fc's that do not reach macroscopic val- 
ues.) By substituting (|138|) in (|137|) and taking the inte- 
gral over the coordinates of all particles but the first one 
we have: 



/oo 
E 
,A fc=l 



{-l)\ kk d k p 



kl 



dz h 



dri, (139) 



where we resorted to Eq. (j28[) to reconstitute the first 
term of the sum. Summing the Taylor series yields: 

lnS^ - InSy = p J [P{ze-P ui ) - P{z)} dr u (140) 
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or, by turning to the ft-potential 



fty — fty 



[P{ze~ 0ui ) - P(z)] dri (141) 



and, by eliminating the same terms, we write: 



P( 



)dn 



= - J P[fi- u(r)]dr = - J P(r)dr, 



(142) 



that is the outcome apprehensible from the point of view 
of physics, expressing the additivity of thermodynamic 
quantities of a nonuniform system. Here, we applied the 
well-know expression for the chemical potential in a field: 



H = fi (P) + u(r), 



(143) 



where fio is the chemical potential observed in the ab- 
sence of a field as a function of pressure. 

Let us make use of relation ()142|) to calculate the 
" 'surface tension" ' of a boundary of macroscopic thick- 
ness. Assume that the potential have a solid core and a 
smooth boundary. For the sake of simplicity we also pre- 
sume that it does not turn negative anywhere and that 
it monotonously falls as it distances away from the field 
region. We shall further denote the core volume by V\, 
and that of the transition zone - by Vt ■ Whereas Vt -C V\ , 
from (|142p we may write 



ft = -P(V — Vi — V t ) — A J P{r)dx, 



(144) 



where A stands for the area limiting the field region, P 
- is the pressure in the region where the field is absent, 
L t - is the length of the transition zone, while x - is the 
coordinate directed along the gradient of the field. Com- 
paring (|144p to the standard definition of the ft-potential 
for a two-phase system in equilibrium 

ft = -P{V - Vi - V t ) + oA (145) 

(since the system volume reduced by V\ + Vt), we have 



(146) 



at the same time, the number of surface particles Ns is 
defined as: 



Ns 



da A . d 



dfi 



= A wJ p(r)dx 

L t 



(147) 



A I q-P[v - u(r)]dx = A I g{r)dx. 



So, in this case we have a < 0, and Ns > 0. 

However we can include the transition zone to the vol- 
ume term just as well and, alongside (|144p . come up with 
an alternative outcome: 



ft 



-P(V -Vi)+A / [P - P{r)} dx. 



(148) 



Provided that: 



a = J [P ~ P( r )} d x 

Lt 



(149) 



and 



Ns = - 



A J[q- g(r)]dx, 

Lt 



(150) 



where g - is the density in the region where the field is 
absent. We can see that now it is a > 0, and Ns < 0. 

Evidently, these two points of view have equal rights in 
this case and are accounted for by two different ways of 
taking surface effects into consideration: either by adding 
them as external ones, or subtracting them as internal 
ones. 

We will further see that when we deal with the po- 
tential of a hard solid, like in the case of the partition 
function of the OSE described above, there are no such 
arbitrary options. 

It is interesting that an attempt to set Ns = by 
shifting the boundary between the volume and surface 
components, like in the case of the phase equilibrium, 
actually proves to be a bad idea. The point is, relations 
(|144[) - (1148)) are true for arbitrary z and T, since they are 
not bound by the phase equilibrium curve. The surface 
obtained on the condition of TVs = will, strictly speak- 
ing, depend on the value of the chemical potential at a 
given temperature and thus loose the physical meaning. 

It is possible, though, to set to zero the coefficient at 
the first degree of activity (density). As it can easily be 
obtained from (|142p this condition is set by 



(151) 



(1 — (pi)dxi = / tpidxi 



where the integration in the left-hand side of the equation 
is performed from the field region to the outside, and in 
the right-hand one - inward the environment, defining the 
location of surface Iq that is of interest to us. At that, 
for the surface tension coefficient and number of surface 
particles the series expansion in powers of activity starts 
from the second degree. 

So, we have three options of considering of surface ef- 
fects: one of them starts with a quadratic term of activity, 
two - with a linear ones. We will discuss them in section 
IVIDI to provide an example, that is closer to reality. 
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C. Sharply Varying Field 



Let us now return to the equation (|137[) and move from 
the smooth field to the opposite extreme case - the field 
of a hard solid. The potential energy of particles may be 
written 



+00 (n e v) 

(r 4 i v). 



(152) 



At the same time, functions ipi (|128p take the form of 
window functions (JTJ) : 



9 1 



(153) 



Let us assume that this body occupies the macroscopic 
volume v inside V . As it has already been stated, we pre- 
suppose that v <C V and the body is far enough from the 
boundaries of volume V. Having set so, we can elimi- 
nate the sign V in the integrals over the coordinates of 
particles from (|137|) and consider the integration to be 
performed over all space: 



In Si 



(-l) fc 



fc=i 



fc! 



- r l...fc a ' !■ 



■ dri 



(154) 

So far, this equation does not allow us to perform inte- 
gration according to (|116[) due to a distortion occurring 
at the boundary of volume v. For example, when the first 
particle is at the distance of a "'correlation length'" from 
the boundary then, in order for (|116[l to fulfill, the inte- 
gration over the coordinates of remaining particles shall 
be carried out over all space, and shall not be limited, 
like in this case, by functions ipf. In order to employ 
formula (|116[) let us reconstitute the integrals in (|154[) in 
the way similar to (f6"6f : 

In - lnS y = £ / ^^dn-dr* (155) 



E 

fe=2 



(-1) 



fc=l 
fc 



fc! 



.1=2 



1 ^fe^i-drfe. 



This treatment will ultimately be analogous to the one 
given above for the OSE, however, we will carry it out 
at the level of correlation functions instead of that of 
factors. 

By taking into account (| 1 16|) . the first sum yields 
f3v[P(0) — P(z)], and in the second one we may express 
i>\ for i=2,3... through \\ &■ 

We will thus obtain: 



lnS^-lnSy = Pv[P(0) - P(z) 



+ E 



k=2 



(-1)' 

fc! 



na-x?)- 1 



.1=2 



(156) 
3 r i k ]k dr i-~ dr k- 



In accordance with our supposition of the analyticity of 
P(z), we assume that P(0) = 0. Once again, as it can be 



easily seen, the structure of integrands in the right-hand 
side of (|156p is such that at least one pair of particles al- 
ways lies at the different sides of the field boundary. At 
the same time, due to the local properties of f[ k \ the in- 
tegrals are defined, like in (|70|) . by the region close to the 
surface. We may once again take the integral along the 
boundary surface over the corresponding coordinates of 
the first particle, the integrands being independent from 
them in this case. Upon making these transformations 
and eliminating the same terms, we write: 



nU = -{V-v)P(z,T)+aa(z,T), 



(157) 



where a - is the area of the surface of a hard solid sub- 
merged in a fluid, 

OO / \ fc n 

P{z,T) = gk B T~ k B T^2 prp- / ^^...dr*, 

(158) 



k=2 



and 



a(z, T) = fc^Ty^ 



00 / 1\fe 



k=2 



(-1)' 

fc! 



(159) 



x / i>i 



1 -f[(i-x0 



1=2 



T ( ^ ) k dxidr 2 ...dr k , 



where x\ - is the coordinate of the first particle perpen- 
dicular to the boundary surface. 

Two equations (|158p . (|159|) are equivalent to the pair 



(|68|) . (|74p in the cases when the latter ones are true. This 
can be easily proved by substituting expansion (|108[) to 
(|158p , (|159p and changing the order of summing. 

The expressions analogous to (|158p . (|159p for a multi- 
component environment may prove effective in the case 
of employing the Coulomb interaction of particles. 

Relation (|157p has the same form as it does for a two- 
phase system in equilibrium, with the exception of the 
fact that cr - the surface tension coefficient - started de- 
pending on two variables. Indeed, upon introducing the 
field of an infinitely positive potential, the system volume 
became equal to V — v and we obtain (|157p . 

Note that, as it has been previously said, we cannot 
move the volume boundary in the case of the field of a 
hard solid, unlike that of a " 'thick" ' surface. The bound- 
ary is rigidly set, whereas the volume and surface terms 
have different functional relationships with the activity: 
the linear relation is opposed to the square one at low 
densities. Thus, in this case, there is no arbitrariness in 
the sign and value of a and N$. 

So, the expression for <j{z, T), occurring in the expres- 
sion for the partition function of the OSE corresponds 
ideally to the surface tension observed at the boundary 
between a hard solid body and a fluid. This correspon- 
dence is by no means coincident and will be later dis- 
cussed. 
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D. Actual Solid 



Let us once more go back to Eq. (|137[) and examine it 
for the case of the potential of a real solid body. Obvi- 
ously, we will consider it constant, so that the situation 
corresponds to the metastable state, once again allow- 
ing us to define the surface tension as a function of two 
variables. 

We will first discuss the model case of the field inde- 
pendent from displacements along the surface of the body 
and then we will generalize it. 

Let us separate (|137|) into the volume and surface 
terms like in (|155p . We shall obtain: 



il = -PV + P J tpidn + a u A, 
v 

where we introduced the notation 

" (-1)' 



(160) 



a u (z,T) = k B T^2 



k=2 



k\ 



(161) 



X / Vl 



1-11(1-00 



1=2 



J : [ k \dxidr2---drk 



where x\ - is again the coordinate directed along the 
gradient of the field and 9 is defined as 



= l-ipi = exp(-/3u;)- 



(162) 



Eq. (|16ip . in the same way to the cases investigated 
above, may be reconstituted by employing (|108[) . Then 
we have: 



oo u 
Z 



a u {z,T) = k B T^2 — 



k=2 



(163) 



1=2 



u[ k ] k dxidr 2 ...dr k . 



The first terms of (1 1631) were calculated in a series of 
works by the authors of [12| on the basis of the topological 
approach from Q , though extending it beyond the limits 
of the potential of a hard solid body. 

Let us reconstitute (|160[) to 

x> 

Q = -P(V -V) - PA J (l-<px)dxx (164) 

— oc 

oo 

+PA / ipidxi +a u A, 



where x' is some arbitrary point located inside the tran- 
sition layer or in its vicinity, that defines volume V . 

We can observe the existence of a certain degree of ar- 
bitrariness in distinguishing between volume and surface 
terms. As a matter of fact, this arbitrariness can only 



exist within the range of the transition layer, i.e. x' is 
situated inside it. It is easy to see, indeed, that in the op- 
posite case the issue is but a banal compensation of two 
identical terms of different signs. Thus, like in the case 
with a " 'thick" ' surface, we obtain at least three options 
for taking surface effects into account: 



Lt 



P{V-Vi-V t )-PA / (l-<p 1 )dx 1 + a u A, (165) 



(166) 



CI = -P(V - Vi) + PA / tptdxi + <t u A, 



n = -P(V-V a )+a u A, 



(167) 



where the denotations are the same as in (|144p , and Vq - is 
the body volume delimited by the surface set by condition 

(HSU). 

Option (|165[) attributes surface effects to a solid body, 
whereas (|166l) - ascribes them to a fluid. Option (|167p 
sets linear surface terms to zero by making use of prop- 
erty (fT5T|) . 

Generally speaking, option (|165[) seems to be prefer- 
able when real adsorption is observed, since it is in agree- 
ment with the ordinary notion of particles adsorbed on 
the surface of a solid body. As we can see, in the first 
order term of z it always shows a positive value of the 
surface particles number. 

Option (|167[) is not suited for the case of a real adsorp- 
tion, because it attributes the molecules adsorbed on the 
surface to the volume of a fluid, being in this case an 
apparently artificial approach. Indeed, the presence of 
a strongly expressed negative potential in the vicinity of 
the surface, as can be easily seen from (II 5 1 [) . shifts the 
position of surface Iq inward the body. 

This option is helpful when the surface potential is 
close to the hard solid one, or, in speaking in general 
terms, when adsorption is absent. 

Option (|166|) is the most exotic one; however, it can 
possibly also be efficient in the case of potentials of un- 
usual form. 

Like in the case of a "thick" surface, the option with 
the condition N s — is inappropriate - for the same 
reasons. 

In the case, when we are able to employ (|138p from 
(|16ip we obtain 



cf u (z,T) 



[e 1 p(z)-p(e 1 z)]dx 1 , 



(168) 



which naturally leads, when substituting it in (|165[1 . 
(fTrJfJj) . to (JI31JI, pggj) . Obviously, the limits of integration 
in p68p can be restricted to the transition zone. 

Now let us consider the situation where the field de- 
pends on displacement along the surface. We shall as- 
sume that the body is a crystalline solid and thus have 
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the surface structure recurring at some periodical inter- 
vals. 

Examining the integrals over the coordinates of the 
first particle that are parallel to the surface obviously 
yields their averaging over such a surface unit cell. Its 
linear dimensions can significantly supersede those of a 
volume cell at certain structural orientations of the crys- 
tal in relation to the surface. 

So, ([T64]) takes the form 

x 

ft = -P(V - V) -PA J (1 - <px)dxi (169) 



PA I (pidxi + a a A, 



fluctuation: 



where (p% - is the value averaged over the surface cell at 
a set depth of a transition zone, and for a a we have: 



Ou(z,T) = k B T^2 



k=2 



(-1)* 

k\ 



(170) 



l-f[(l -ft) 



i=i 



f\ .k dr 2--- dr k) 



cell 



where the triangular brackets also stand for the averaging 
over the surface cell at a set depth of a transition zone. 
We should remind our reader that !F\ is not directly 
subject to averaging, as they belong to the unperturbed 
environment. 

The counterpart for (|163p will be 



k=2 



Z 



(171) 



dxi( 0i 



l-ii(i-w) 



1=2 



U[ k) k dr 2 ...dr k ) 



cell 



where, once again, averaging is performed over the sur- 
face cell at a set depth. 

Equations (| 1 5 1 [) . f| 1 65|) - (|167p maintain their forms, 
with ipi — > ipi and a u — > a a being substituted. 

In the case, where the contact surface contains crystal 
faces differently oriented in relation to the crystal axes, 
the surface terms of expression (| 169[) should be summed 
over these faces. 



Po oc exp {-(3R min ), 



(172) 



[lol , p. 339], where R m in - is the minimal work done on 
eliminating the fluctuation. 

Indeed, as the chemical potential for the particles 
located outside of the field does not change upon a 
local introduction of a field according to (|143[) . then 
we may assume that, in this case, the process of cre- 
ation/elimination of a fluctuation is carried out at a con- 
stant chemical potential and temperature, i.e. the natu- 
ral variables of the Sl-potential. Thus, the minimal work 
on forming a cavity is equal to the alteration of the f2- 
potential: 



Pi, 



= n 11 - n = [-P(v - v) 



-PV] 



= vP(z,T) + aa(z,T). 



(173) 



As we are interested in the fluctuation that causes den- 
sity to change in a stepwise way, this process of creating 
the fluctuation must be carried out by a hard solid. It is 
here where we find the meaning of the correspondence ex- 
isting between the surface tension coefficient in the OSE 
and that for the boundary surface of a hard solid and a 
fluid. 



B. The Difference between the OSE and the GCE 

The core difference between the two ensembles is in 
surface effects, that can be observed in (|63|) . In the case 
where the system does not interact with the environment, 
Sy = E] v E]v-v and this equation is deduced to a banal 
one 

lnT„= lim [In (E v Ev- v ) - 1nE V -v] = lnH„, (174) 

V^oo 

demonstrating the identity of the ensembles. 
At v = V from (16111 we obtain 



InPn 1 



In 



(175) 



just like it should be from the point of view of an ordinary 
treatment of the GCE. 

However, the differences are not connected to the pre- 
sumed lack of surface terms in the GCE; they are due 
to fact that such terms have another form. Indeed, by 
applying the operation of singling out the surface term 
to (162[) we may write: 



VII. DISCUSSION 



The Implication of the Surface Term for the 
OSE 



Eq. (|75p for the probability of cavity formation, which 
defines the partition function of the OSE, is in agree- 
ment with the general expression for the probability of a 



In Si 



E 

oo 



t=2 



V 



X 



V\ 



(176) 



U± ' t dr\...dr t . 



By comparing (I176|) with (l67|) . (|73|) we witness that 
Z v = expj3[VP(z,T)-Aa(z,T)], (177) 
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where we resorted to the substitution tp <-> \ m the last 
term of (|176p , which always remains possible due to the 
spatial symmetry of the task. 

Let us write the formulas (72]), dUJ, (flUO)) of the OSE 
once again for the sake of comparison: 



T v =expP[VP(z,T)+Aa{z,T)}, 



dlnz 



N b -N S =N- N 8 



and, by differentiating (jl77[) . those of the GCE: 



d In E v 
dlnz 



= N b + N s = N. 



(178) 



(179) 



(180) 



(The last equality in (|180[) . as is known, can be easily 
obtained by differentiating (p~0|) over the chemical poten- 
tial.) 

Thus, the distributions of the OSE and GCE contain 
the surface terms of different signs. We will not be ex- 
amining this issue profoundly right now. We would like 
to briefly note, though, that the reverse sign of the sur- 
face term in (I177P is no coincidence, and that the GCE, 
being a set of closed systems, cannot reproduce the open 
system in all its completeness. 

Now we can turn our attention back to (f2"9"| . From (TTF 
and ([M)]) follows that 



1% Q%\ty V )dri = N = N b + N s 



(181) 



therefore, we see that Qq\ k contains surface terms and 
that is the exact reason for (|29p to not be fulfilled when 
the first particle is located close to the boundary of vol- 
ume V. 

As a comparison, from (|126p and (|179p we have for the 
OSE: 



^(fpdrx =N = N bl 



(182) 



which is due to the exact satisfaction of (| 1 1 6[) indepen- 
dently of the position of the first particle. 

Eq. (|175p demonstrates the lack of scalability of the 
GCE or, in other words, it shows that this distribution 
does not match the conformity criterion. 

The GCE does not scale because in the case of v = V it 
yields the expression different from the case where v < V, 
which is evident from the comparison of (|6ip and (|175p . 
We shall further note that the volume term, obviously, 
maintains its form and, thus, this part of the GCE is 
scalable. 

Indeed, considering expression (|58p . leads to the un- 
derstanding that the presence of the window functions 
ipY distorts the correlation functions for those configu- 
rations when a group of particles from the 1st to kth is 
located in the vicinity of the boundary of volume V. It 
happens as a consequence of the integrals of the Ursell 
factors do not reach in this case their extreme values as in 



integrating over the whole space or the whole area where 
they are different from zero. 

Clearly, the correlation functions in the uniform envi- 
ronment should possess the property of translation in- 
variance, which is not fulfilled for GCE in the vicinity of 
the boundaries of volume V. In order to eliminate this 
distortion we should pass to limit (|89p . and that is ex- 
actly what happens particularly in equations (|63p . (|90p . 
It is by liquidating this distortion that we can obtain the 
right expressions for the surface tension in the OSE. 

Let us also note, that such a transition is impossible if 
we make use of relations like (I15p . incorporating extensive 
values. However, they can always be deduced to form 
(f24|) . by, for example, employing (|20p. and after that we 
turn to limit 



C. Small Fluctuations 

Logical reasoning tells us that the value of surface 
terms should decrease for the terms of the distribution 
of the OSE approaching the average quantities. It would 
be strange indeed if the probability to find an average 
number of particles in a certain volume contained the 
surface. As it has already been said such compensation 
does take place. 

First of all, the very equation (|100p indicates that. 
Since, at small values of the particle number in a given 
volume, the probability depends on the size of the limit- 
ing surface, which is clear from (|75p . and since the aver- 
age value does not depend on it, the only option would 
be compensating the surface values for the terms of the 
distributions with larger m's. 

Second, according to the calculations, we observe the 
strict compensation of the surface terms in the vicinity 
of average values for the first terms of the expansion over 
activity. 



VIII. SUMMARY 

The article consists of two parts: the description a 
new statistical ensemble and the consideration of surface 
effects for a system in the field that corresponds to a solid 
body. Although, the latter part is auxiliary it contains 
some results that are of their own individual value. 



New Ensemble 

1. We introduce a new statistical ensemble named the 
OSE (open statistical ensemble), whose major fea- 
ture is in correctly taking into account the surface 
terms for an open system, which leads to substitut- 
ing the Boltzmann factor in the configuration inte- 
grals of the grand canonical ensemble (GCE) with 
the correlation function of a specific form (|118p . 
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2. The partition function of the OSE in its expanded 
form is defined by expression (|64p , which includes 
volume and surface terms (|72| . The analogous 
equation for the distribution of the OSE is given 
by the sum ([90]) , whose first term represents GCE. 

3. Unlike in the case of a two-phase system in equi- 
librium, the surface tension coefficient for the OSE 
depends on two variables: pressure as well as tem- 
perature (f?4"]) . 

4. The expression for the partition function of the 
OSE conforms with the thermodynamic approach 

una). 

5. Thermodynamic and statistical relations of the 
OSE (sectionfVj are ensured by the recurrence rela- 
tions (|5ip - (156)) of a new class of correlation func- 
tions B™ an d by the use of generating func- 
tions (subsection IV Dp . 

6. Functions (|4l))) generalize the notions of 
the Boltzmann (|41[) and Ursell (|4"2"]) factors and con- 
tain them as extreme cases (|4"8"j) . (|49l) . 

7. Functions B[ m ' 7 ^ +k and A[ m '^ +k contiguous to 
them are yielded by fractional generating functions 

(®, dSHD- 

8. The correct consideration of the integrals of Ursell 
functions for the OSE results in obtaining the exact 
analogs for the Ornstcin-Zernike equation of higher 
orders (|116[) . (| 1 1 2[) is a further generalization of 
this relation. 

9. The basic equations of the OSE can be obtained 
on the level of factors as well as on the level of 
correlation functions thus alleviating the problems 
related to the convergence of series in density. For 
instance, we can thus yield the expressions for the 
surface terms (|161l) and (|163|) . 

10. The mean particle number for the OSE does not 
contain surface terms (llOOp . which indicates them 
being compensated by the high-order terms of dis- 
tribution. 

11. The GCE contains the surface terms, whose sign is 
contrary to that of the OSE (flTT]) . ([178"]) . 

12. As a consequence, many expressions of the GCE 
feature inaccuracies related to the surface effects 
HMD, (gEED - 



A System in the Field 

1. The investigation of the system of particles in a 
field of force of general type proves that the surface 
terms of the OSE (|159p really correspond to the 
surface tension observed at the boundary separat- 
ing a hard solid from a fluid. This rigid character 
of the solid is explained by the fact that from the 
viewpoint of probability distribution we are inter- 
ested in the configurations abruptly breaking off in 
density. 

2. The consideration of a solid body as an external 
field is a productive approach, serving as an alter- 
native to the equilibrium one (a nonvaporizing or 
insoluble body). 

3. The statistical approach gives a closed expression 
for the fl - potential of a nonuniform system (|142p 
in the case of the macroinhomogeneity, which has 
a clear physical meaning. 

4. In general, the volume term, adsorption and surface 
tension cannot be separated with outmost precision 
(|160p . and it is expedient to interpret them on the 
basis of the form of the potential of the body and 
fluid particles interaction (|165p - (|167p . 

5. In particular, depending on the situation, surface 
terms can start with both a linear and quadratic 
term of activity (US]), dUTJ. 

6. The condition N s = is, by all appearances, al- 
ways adverse for the case being considered, since 
the position of the separating surface (]15ip starts 
depending on pressure. 

7. Surface terms, like in the case of a sharply varying 
potential, can be calculated in its general form for a 
smooth field (|168p as well as for a real one (I163P , in- 
cluding the case of the surface inhomogeneity (|171[) . 

8. In the case of a smooth field the series for the sur- 
face term converges to a compact form (|168p . 
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